Interference scheme to measure light-induced nonlinearities in Bose-Einstein condensates 
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Light-induced nonlinear terms in the Gross-Pitaevskii equation arise from the stimulated coherent exchange 
of photons between two atoms. For atoms in an optical dipole trap this effect depends on the spatial profile of 
the trapping laser beam. Two different laser beams can induce the same trapping potential but very different 
nonlinearities. We propose a scheme to measure light-induced nonlinearities which is based on this observation. 
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In the framework of the mean-field theory an atomic Bose- 
Einstein condensate (BEC) represents a nonlinear system 
where the nonlinearity is caused by atomic collisions. Vir- 
tual photon exchange in a BEC of neutral polarizable atoms, 
caused by off-resonant optical laser radiation with the Rabi 
frequency and the detuning A leads to additional nonlin- 
earities in the condensate's equation of motion d 0, HI H Si • 
Since the optically induced nonlinearity (OINL) is strongly 
suppressed for large detuning, it can be safely neglected for 
atoms trapped in optical dipole potentials |6]. However, for 
smaller detunings it must be taken into account. It has been 
theoretically shown that the OINL in a BEC can lead to 
nonlinear diffraction of ultracold atomic beams laser- 
induced gravity |5], and one can create Thirring and gap soli- 
tons 1 1], photonic band gaps and defect states in a periodic 
BECH. 

Inspite of the fact that the OINL can be of the same or- 
der of magnitude as collisional nonlinearities |4], the obser- 
vation of the OINL is difficult. Due to spontaneous emission 
an atomic BEC in a largely detuned laser field typically deco- 
heres at a rate of 7d ec := 7|^| 2 /A 2 , where 7 is the sponta- 
neous emission rate of a single atom in free space (see, for in- 
stance, Ref. 1 8]). Hence, in order to avoid spontaneous emis- 
sion one has to consider only short interaction times T such 
that T7dec < 1. In this case, the effect of the OINL essen- 
tially amounts to a small nonlinear phase imprint on the BEC. 
The goal of this paper is to develop a scheme to measure this 
phase imprint. 

We would like to note that decoherence of another type can 
be caused by the fact that there is a distance (the so-called 
Condon distance) at which a pair of atoms is resonant with 
the light, even if the light is detuned from the atomic reso- 
nance for a single atom |9]. However, we neglect this type of 
decoherence, because in the case of 87 Rb at the typical density 
of 10 14 cm -3 and the detuning A ~ 1 GHz, the Condon dis- 
tance is almost five times smaller than the mean interatomic 
distance. 

The detailed derivation of the OINL in a BEC was given in 
Ref. 0]. Its features are determined by the properties of the 
local electric field, which governs the evolution of the con- 
densate's wave function ip in the field of electromagnetic ra- 
diation. In the case of low light intensity and in the slowly 
varying (in time) amplitude approximation, we have the fol- 



lowing equation for the local electric-field strength 0,Q3l 

E+ C (x) = E+ (x) + J dx'e^«P+(x') 



(1) 



3 cos 2 6-1 ,3 cos 2 6»-l 



R 3 



R 2 



cos 2 6-1 
R 



where R = |x - x'|, P+(x) = a\ip(x)\ E ; + C (x) the 
medium polarization, a = — d 2 /(HA) the atomic polariz- 
ability, and 6 the angle between the vector R and the dipole 
axis. in Eq.Q is an incident laser field with the fre- 

quency ojl = ckr,, which is close to the frequency ljq of the 
electric-dipole transition. The integral term describes dynam- 
ical dipole-dipole interaction and leads to the OINL in the dy- 
namics of the condensate's wave function ip. As it was shown 
in Ref. |4], Eq.Q describes subsequent pair interactions be- 
tween the atoms when an off -resonant laser photon enters the 
medium of ground-state atoms and virtually excites an atom, 
then this atom goes back to the ground state and emits a pho- 
ton, which is absorbed by another atom and so on. We would 
like to stress that the mechanism described above is different 
from the (quasi)photoassociation, when a pair of atoms first 
absorbs a photon and undergoes a virtual transition to an elec- 
tronically excited quasimolecular state and then it reemits the 
photon and returns to the initial electronic state fllll . 

The dipole field in Eq.Q contains three terms, which have 
different dependence on R. The first one (^ 1/ i? 3 ) is a near 
field and gives a significant contribution at distances much 
smaller than the laser radiation wavelength. The last one 
(~ 1/R) is a far field and gives a dominant contribution at 
distances greater than the wavelength. At distances of the or- 
der of wavelength all the three terms are of the same order of 
magnitude. They must be taken into account for condensates 
with a typical size of several optical wavelengths. 

The near and far fields have different dependencies on the 
angle 6. The sign of the near-field term can be either negative 
or positive depending on 6 and the near field vanishes after 
the angular integration in the case of a homogeneous density 
distribution. The sign of the far field does not depend on 6 
and it does not vanish for any atomic distribution. This leads 
to an attractive effective nonlinear potential in a large system 
of interacting optically induced dipoles. 

The above discussion refers to the microscopic electrody- 
namics of a BEC. Since usually the condensate density does 
not change significantly over one wavelength, the nonlinear 
optical potential can also be described in terms of a macro- 
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Our scheme to measure the OINL is based on the obser- 
vation that K, pt ,±(x) can induce about the same mechani- 
cal forces on atoms even for very different laser beam pro- 
files. A doughnut laser beam with Rabi frequency rido(x) = 
f2° exp(— r\/w 2 ) exp(ikLz)(x + iy)/w and a Gaussian laser 
beam with ^Gauss(x) = f2 Gauss exp(-r^/w 2 ) exp(ikLz) pro- 
duce the optical potentials 
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FIG. 1: Level scheme of an atomic spin-1 BEC interacting with 
off-resonant circularly polarized laser waves. 

scopic (averaged) electric field |4]. This approach is equiva- 
lent to the microscopic one and allows to simplify the equation 
for the condensate wave function. As a result of the macro- 
scopic approach the OINL can be written as a local attractive 
term [see Eq.@] in the Gross-Pitaevskii equation (GPE). 

In this paper, we are concerned with atoms possessing two 
internal ground states characterized by the magnetic quantum 
number m = ±1. In the presence of two circularly polarized 
laser beams and in the case of a dilute gas with the level con- 
figuration shown in Fig. ^ the two collective wave functions 
V>± fulfill the modified GPE fl 
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+ [K s + K a + K pt,±(x)] \lp±\ + (k s - K a ) |<0 T | >1p± , 

where M is the atomic mass. It is assumed that the detunings 
A± of the laser frequencies from the respective atomic tran- 
sitions satisfy the inequality |A_| ^> |A+|, such that there is 
no population exchange through Raman transitions between 
the components and the condensate is similar to a two-species 
one 0| . The strength of the collision-induced nonlinearity 
is described by the coefficients n SM := Airh 2 a s ^ a /M, where 
a s<a are symmetric and antisymmetric scattering lengths, re- 
spectively. For large detunings (|A±| 3> 7) one can neglect 
the refractive index of the gas even for interatomic distances 
of the order of the transition wavelength Ao. In such situa- 
tion the Rabi frequencies f2±(x), which are proportional to 
the electric-field strength of the a T -polarized laser beams, do 
not depend on the atomic density and the potential 
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does not lead to any nonlinearity in Eq.(|2|i. The OINL coeffi- 
cient is also, in general, position dependent and in the case of 
low light intensity (|f2±| 2 <C A±) is given by 



«o P t,±(x) = -27T?17 
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r± := \fx 2 

the same harmonic trapping potential with trap frequency u> 
provided that 



y 2 . Thus, apart from a constant energy shift 
2 /Acauss, both beams provide approximately 
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holds. Note that Ado must be positive and Acauss negative in 
order to provide a trapping potential. 

A very important physical difference between these two po- 
tentials is that in the doughnut potential atoms are low-field 
seekers, whereas in the Gaussian beam they tend to accumu- 
late at positions with high laser intensity. This implies that the 
atoms experience a strong OINL in a Gaussian beam, but a 
negligible nonlinearity in a doughnut beam. This effect is the 
fundamental phenomenon on which our proposal is based on. 

We would like to emphasize that in this situation any differ- 
ence in the evolution of the two BEC components |±) would 
be produced by the OINL alone. In the harmonic approxi- 
mation, the trapping potentials only differ by an unimportant 
spatially homogeneous value Sy- Furthermore, for 87 Rb the 
antisymmetric scattering length a a is negligible (see numeri- 
cal values below), so that the collision-induced nonlinearities 
are identical for both components. 

The measurement scheme consists of the following steps. 

(i) 7r/2 Raman pulse preparing a superposition of the states 
|+) and \— ). We assume that the BEC is initially prepared 
in the state |— ) with a spatial wave function ^>(x, 0) corre- 
sponding to the stationary collective ground state in a dough- 
nut beam with er + polarization (fi_ = fido)- The detuning of 
the doughnut beam is huge so that both the spontaneous emis- 
sion and the OINL can be neglected. Half of the atoms are 
then transferred to the state |+) by a 7r/2 laser pulse. Since 
powerful experimental techniques such as stimulated Raman 
adiabatic passage 11211 do exist, we can assume that the trans- 
fer is nearly perfect and can be described by the unitary oper- 
ator 
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where ip s is the phase of the <j_ -polarized Stokes laser. It can 
be freely chosen and will be discussed later. After this trans- 
formation with the duration TL/2, the BEC is in a coherent 
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superposition of the states |±) with the wave functions 



N-/N 



v/2 
V>(x,0) 
v/2 



(9) 



(ii) Phase imprint. In the second step, we apply an addi- 
tional <r_ polarized Gaussian laser beam to the BEC so that 
fl+ = r^Gauss- Both detuning and intensity of this beam are 
much smaller than that of the doughnut beam and are cho- 
sen so that relation Q holds. For small enough detuning the 
OINL is not negligible for this beam but completely negligi- 
ble for the doughnut beam. Thus, only the wavefunction for 
atoms in state |+) is changed by the OINL. Since all other 
parameters are held constant this is the only change the sys- 
tem experiences beside the phase shift caused by the potential 
difference Sy- 

Since the interaction time T must be small as compared to 
l/7dec it has to be of the order of 10 /is in a realistic experi- 
ment. For such a short interaction time, the spatial distribution 
of rubidium atoms does not change much and, therefore, it is 
possible to neglect the center-of-mass motion. The wave func- 
tions at time t = T + then can be written in the form 



(10) 



V2 



where (^oinl = — K pt,+(x) \ip+\ T/Ti is the phase shift pro- 
duced by the OINL, i.e., the quantity we are interested in. /j, 
is the chemical potential and cpy ■= SyT/h. 

(iii) Second tt/2 Raman pulse. After the phase imprint 
the Gaussian laser beam is switched off. The interferometric 
scheme is completed by applying a second tt/2 Raman pulse 
that is identical to the first one. After this pulse a number of 
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atoms are left in the state |— ), where they are trapped by the 
doughnut beam. The atoms in the state |+) are not trapped 
anymore and are dispersed in a time that typically is of the or- 
der of a few milliseconds. To make N- a signal of the OINL 
alone, one has to eliminate the influence of the potential dif- 
ference Sy which can be done by choosing ip s = —ipy/2. In 
this case, N- vanishes in the absence of the OINL. 

In order to get concrete estimations, we have considered 
a BEC of N = 10 5 rubidium atoms ( 87 Rb, M = 1.45 x 
1(T 25 kg, a s = 5.4 nm, a a = -0.05 nm 0, A = 780 
nm, 7 = 38 MHz) in a trapping potential of the form (|5) 
and with the parameters A do = 1.1 x 10 15 Hz, fig = 
3.15 x 10 10 Hz (which corresponds to the oscillation fre- 
quency ujj_ = 2ir x 576 Hz), w — 10 nm. Using these 
parameters, we have performed analytical two-dimensional 
(2D) calculations along the lines described above by em- 
ploying the Thomas-Fermi approximation for the calculation 
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FIG. 2: Number of atoms in the trapped state | — ) after the inter- 
action with the Raman pulses. Dots - full numerical simulation, 
solid line - numerical integration of Eq. Jl H with exact ground state, 
dashed line - analytical calculations based on the Thomas-Fermi ap- 
proximation rEq. J12H . 



of the ground state wave function ^(x, 0). The condensate 
has been assumed to be homogeneously extended along the 
z direction with a length of L z — 20 /_tm, which gives 
the following value of the Thomas-Fermi radius: Rtf — 
2 (h/Muj ± ) 1/2 [(a s + a a )N/L z } 1/4 = 2 fim. The number 
of atoms remaining in the trapped state iV_ can be worked 
out according to Eg. (II 1> . assuming that the condensate size is 
much smaller than the laser width w. The result is given by 
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and is shown in Fig. |2] (dashed line) for the interaction time 
T = 10 /is. As one clearly sees, the number iV_ of trapped 
atoms increases with the increase of the light-induced nonlin- 



earity which is proportional to 
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up to 0.25N could 



fill T7dec <§C 1 the parameter \Qq. 
higher than about 0.001, i.e., values of N. 
be observable. 

To confirm our analytical results, we have also performed 
2D numerical calculations in which the center-of-mass mo- 
tion, the full form of the potentials Kpt.±, and all nonlinear- 
ities are taken into account. The initial ground-state wave 
function ip has been obtained by numerical solution of the 
two-dimensional Gross-Pitaevskii equation. For the time evo- 
lution, we have employed the split-step method (see, e.g., 
Ref. 11511 ') and combined it with the imaginary time propa- 
gation technique (see, e.g., Ref. flrjlo to find the ground state. 

Starting from state (|9jl, we have solved the time evolution 
of the BEC which is governed by Eq.©. After that the trans- 
formation U of the two-component wave function has been 
performed and the number of atoms in the trapped state A r _ 
has been worked out. The dependence of N- on the param- 
eter |^Gauss| /^Gauss ^ s shown in Fig|2] (dots). The results 
are in good agreement with the approximate analytical calcu- 
lations. However, the Thomas-Fermi approximation leads to 
somewhat higher values of compared to the exact numeri- 
cal calculations, because in neglecting the center-of-mass mo- 
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tion the Thomas-Fermi approximation gives a slightly higher 
central density. The comparison of our analytical and numer- 
ical calculations shows that the center-of-mass motion does 
not influence significantly iV_, which is mainly determined 
by the form of the initial atomic distribution and the param- 
eters of the laser radiation in the middle of the trap. After 
the second Raman pulse those atoms transferred to state |+) 
are dispersed within about 1 ms, because the Gaussian laser 
is switched off. Since the transfer is not adiabatic the num- 
ber of atoms in the state |— ) is reduced, but the shape of the 
wave function remains almost the same. If there would be 
no atomic collisions this would be again a stationary state of 
the corresponding Schrodinger equation, but in the presence 



of collisions this is not the case and the remaining cloud of 
atoms in state |— ) starts to oscillate radially. They can be de- 
tected by taking an absorbtion image with resonant laser light. 

In the present paper, we have developed an interference 
scheme that allows to measure the OINL caused by dynam- 
ical dipole-dipole interaction in the field of an optical laser. 
The observation of this effect would help to improve our un- 
derstanding of the interactions between light and atoms. 
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